﻿ Transformata Laplace Transformata Laplace generalizează ideea transformatei Fourier in tot planul complex Pt un semnal x(t) spectrul sau transformata Fourier este ∞ − jtω Xxte dtω= () () ∫ −∞ Pt acelaşi semnal x(t) se poate introduce transformata Laplace ∞ st− () ( ) ;Lxt s X s xte dt== sjσω=+ () () () ∫ −∞ ∞∞ = de la − la unde integrala se efectuează pe o curba cstσ = transformata Laplace este spectrul semnalului adică transformata Fourier Pentru 0σ jω σ=cst σ Domeniul de convergenta = domeniul valorilor lui s pt care integrala este convergenta sau X(s) - finita Proprietati DC 1 format din benzi paralele cu axa imaginara jω ine nici un pol al transformatei Laplace 2 nu conţ 3 pt semnale cu întindere spre dreapta (o subcategorie fiind semnalele cauzale), DC începe cu cel mai dreapta pol si se întinde tot spre dreapta { } maxReps >σ 4 pt semnale cu întindere spre stânga (o subcategorie fiind semnalele anticauzale), DC începe cu cel mai dreapta pol si se întinde tot spre stânga { } minReps axa jω inclusa in DC 1 Transformata inversa ∞ 1 st xtXseds= ()() ∫ 2jπ −∞ ) Integrala se efectuează pe dreapta paralela cu axa imaginara (cstσ = Mai des se foloseşte metoda descompunerii in fracţii simple si inversarea cu ajutorul tabelelor Probleme 2s+ 1 Transformata Laplace X(s) a unui semnal x(t) este ()Xs= 2 (1)(3)ss− + Sa se determine x(t) in următoarele condiţii a) Re(s) 1 Rezolvare 12AA2Bs+ ++ ; 22()Xs== (1)(3)s-1 s+3s-1ss−+ () 2s+ = =− => (3)() 1/16Bs Xs=+ 23s=− (1)s− 3s=− 22s+ = = 2(1)() 3/4As Xs=− 1s= + 13ss= − 232dss+− =− = = 12(1)() 1/16AsXs⎡⎤ ⎣⎦ + 11(3)ssds s= = 2 1/16 3 / 4 1/16s+ +− 22()Xs== (1)(3) s-1 s+3s-1ss−+ () 313 1tt t− − − − + −σσ σ ()xtettetet=− a) DC : Re(s) () ()() 16 4 16 313 1tt t− σ=− − − − − b) DC: -3 () ()xtettetetσσ ()() 16 4 16 313 1tt t− c) DC: Re(s) > 1 => () σ=+− ()xtettetetσσ ()() 16 4 16 Verificati expresia lui x(t) folosind mediul Matlab: >> L=(s+2)/((s-1)^2*(s+3)) >> ilaplace(L) 2 2 Funcţia de transfer a unui sistem liniar invariant in timp si cauzal este data de expresia 2s+ ()Hs= Determinaţi si desenaţi răspunsul sistemului, daca la intrarea sa se 2 22ss++ 2t− aplica semnalul ()xte= Rezolvare ytxtyt YsXsHs=∗↔ = ⋅ () () () ( ) ( )() 22ss++ ()Hs== cu poli complecşi 1sj= −±1⇒σ =− 21,2maxp − − 1222( )( )ss ssss++ jω DC pentru H(s) 1 σ -2 -1 -1 Sistemul este cauzal: DC este { Re 1s >− } 2t− ⎧ ≥ 222,0tttet− − == = + − ()xte e te tσσ () ( ) 2t⎨ ,0et ()Xs=− 2 22 4ss s+− − 42 4s−+ − Ys Xs Hs=⋅= = () () () 22 2 422 222sss sss−++ −++ () ()()() Domeniul sau de convergenta este intersectia DC pentru ( Hs si DC pentru()Xs: ) 1Re 2s− () ()( ) => s= 2, A = -2/5 s=0, C = 8/5 s=1, B = 2/5 ⎤ ++ 21 2 4 21 2 1 3ss⎡ => () + =− + + ⎢⎥ 222Ys=− 525 22 52511 11sss sss−+ + − ++++ ⎢⎥ ()() ⎣⎦ 222ttt−− ⎡⎤ () cos 3 sinyte t etettσσ =−+ + => () () ⎣⎦ 55 3 3 Fie un sistem liniar, invariant in timp si cauzal a cărui funcţie de sistem H(s) are constelaţia de poli si de zerouri din figura: jω σ -3 -1 1 2 i in fiecare caz daca a) Sa se indice toate regiunile de convergenta posibile; b) Indicaţ sistemul satisface sau nu condiţiile de stabilitate si/sau cauzalitate Rezolvare Re(s) sistem anticauzal si instabil, -3 sistem instabil, -1 sistem cauzal si stabil 4 Sa se demonstreze următoarele proprietăţi ale transformatei Laplace: a) deplasarea in domeniul timp; b) deplasarea in domeniul s; c) scalarea in timp x(at); d) convoluţia in domeniul timp; e) derivarea in domeniul s 5 Fie un sistem obţinut prin conectarea in cascada a doua sisteme liniare si invariante in ia de transfer a sistemului timp ale căror funcţii de transfer sunt H1(s) si H2(s) Daca funcţ te sistem invers sistemului echivalent H(s) este unitara, sistemul descris de H2(s) se numeş descris de H1(s) a) Sa se determine relaţia de legătura intre H1(s) si H2(s) toare unui b) In figura este prezentata constelaţia de poli si zerouri a lui H1(s), corespunză sistem cauzal si stabil Sa se determine constelaţia de poli si zerouri pentru funcţia de transfer a sistemului invers asociat spunsul la impuls al sistemului invers in ipoteza ca este stabil c) Sa se determine h2(t), ră d) sa se demonstreze ca răspunsul la impuls al sistemului echivalent h(t) este identic cu δ(t) jω σ -3 -1 1 2 6 Un SLIT cu intrarea x(t) si ieşirea y(t) este descris de ecuaţia diferenţiala: 2 dy dy 23yx+−= 2 dt dt 4 i H(s), funcţia de transfer a sistemului Schiţaţi constelaţia de poli si a) Determinaţ zerouri i h(t) daca – sistemul este stabil; - sistemul este cauzal; - sistemul nu b) Determinaţ este nici stabil nici cauzal Rezolvare 2 2dd ;xtsXs xtsXs↔↔ a) () ( )() ( ) 2 dt dt 211AB = => () => 23sY sY Y X+− =+ 2Hs== −−+−+ 23 1 3 1 3ss s s s s+ ()( ) 11 = = (1)()As Hs=− 1s= + 134ss= 11 = =− (3)()Bs Hs=+ 3s=− − 314ss=− 11 1⎛⎞ Hs=− () ⎜⎟ 41 3ss−+ ⎝⎠ jω σ -3 1 311tt− b) sistem stabil DC: -3 DC: Re(s)>1, () ()() 44 Verificati expresia raspunsului la impuls folosind mediul Matlab: >> L=1/4/(s-1)-1/4/(s+3) >> ilaplace(L) 311tt− sistem instabil,anticauzal DC: Re(s) -2, sistem cauzal Hs== = () ++ ++ 1422 24Xs s s s s s+ ()( ) ()()() () 33jω+ () HHsω== () () = ωω 24sjjjω ++ ()() 33sAB+ () Hs== b) () + 24 2 4ss s s++ + + ()() 333s+333s+ ()() 2()As Hs=+4()Bs Hs= += = = =; () () 2s=−4s=− 42s+22s+ 2s=−4s=− ⎛⎞ 2431 1 3tt−− Hs ht e e tσ =+↔=+ ()()() () ⎜⎟ 22 4 2ss++ ⎝⎠ Verificati expresia raspunsului la impuls folosind mediul Matlab: >> L=(3/2)/(s+2)+(3/2)/(s+4) >> ilaplace(L) M k kk kNMbs∑ Ysdy dx () 0k= c) Ecuaţia diferenţiala () == ∑∑ kkNkkabHs=↔ () 00kkkdt dt X s== kas∑ 0k= ′′ ′ ++= + => 6839yyyxx′ d) Se foloseşte una din cele 2 forme canonice de implementare: ′ ′′′ ++ =++ 01 2 01 2ay ay ay bx bx bx′′ ++ = ++ 012012ayayay bxbxbx′′ ∫∫ = ++ 012012ayayaybxbxbx++ ∫∫ ∫ ∫∫ ∫ 1 => ybxbxbxayay⎡⎤ +−− 0120 1=+ ∫∫ ∫ ∫∫ ∫ ⎣⎦ 20a 6 Forma canonică 1 Forma canonică 2 1/a2 y x b2 x b21/a2 y ∫∫∫ b1 -a1 -a1b1 ∫∫∫ b0 -a0 -a0b0 x 1 ∫ y -6 3 ∫ -8 9 8 Pentru sistemul liniar si invariant in timp cauzal cu răspunsul in frecventa 7jω+ = Hω () 2 41jjω ωω +−+ () () a) determinaţi h(t), b) daţi o structura de implementare constând din doua sisteme conectate in cascada, c) daţi o structura de implementare constând din doua sisteme conectate in paralel Rezolvare 7sABsC++ Hs Hω== =+ () ( ) 22 = +++ 4141sjssssssω +++ () () 73s+ 4()As Hs=+ = = () 2 4s=− ++ 4113sss=− s=0 => C=22/13 => B=-3/13 7 ⎛⎞ 22 22⎛⎞ −− ss⎜⎟ 31 31⎜⎟ 33⎜⎟ =− () ⎜⎟ 22Hs=− 13 4 1 13 413sss s⎜⎟ +++ + ⎛⎞ ⎜⎟ s++ ⎜⎟ ⎜⎟ 24⎝⎠ ⎝⎠ ⎝⎠ ⎛⎞ 22 11 + s⎜⎟ + 31 322⎜⎟ + () 22Hs=− 13 413 13s⎜⎟ + ⎛⎞ ⎛⎞ ss++ ++ ⎜⎟ ⎜⎟ ⎜⎟ 24 24⎝⎠ ⎝⎠ ⎝⎠ ⎛⎞ 13⎜ ⎟ s+ ⎟ 3 1 47 2⎜ 22 + () ⎟ 22Hs=− 2⎜ 13 4 6313s+ ⎛⎞ 13⎛⎞ ⎜ ⎟ s⎛⎞ s++ ++ ⎜⎟ ⎜⎟ ⎜ ⎟ 24⎜⎟ 22⎝⎠ ⎝⎠ ⎝⎠ ⎠ ⎝ 4/2 /23 3 3 3 47 3tt t−− − cos sinht e t e t t e t tσσ σ=− + ()()()() 13 13 2 13 233 Verificati expresia raspunsului la impuls folosind mediul Matlab: >> L=(s+7)/(s+4)/(s^2+s+1) >> ilaplace(L) 717ss++ b) () = 22() ()abHs H sH s== + 4141ssssss++ +++ () () +=+ += Pentru Ha(s): sau 4ay ay bx bx y y x′′ 01 01 ′ ′′′ ′′′′ ++ =++ ++=+ Pentru Hb(s): sau 7ay ay ay bx bx bx y y y x x′′ 01 2 01 2′ x 1 u 1 ∫∫ y -1 1 -4 u ∫ -1 7 8 22⎛⎞ s− 31⎜⎟ 3 c) () =− =+ () () ⎜⎟ 2Hs H s H sαβ 13 4 1sss++ + ⎜⎟ ⎝⎠ Hα(s) x y Hβ(s) 9 Fie sistemul din figura: R R I1(s) I3(s) x C y C X(s) Y(s) I2(s) ia de transfer a circuitului, daca R=100kΩ, C=10uF a) sa se determine funcţ b) Ce tip de caracteristica realizează acest circuit? Sa se determine factorul de calitate Q, amplificarea si pulsaţia de rezonanta naturala ω 0 spunsul sistemului la un impuls unitate, respectiv la un semnal c) Sa se determine ră = armonic () 00cosxtA tω Rezolvare Ys () ; () 3()IssCYs== 1 sC 11⎛⎞ ⇒=+=+ ()()() ()( ) ( ) () 23 2311Is I s R I s I s sCR sCR sCY s=+ ⎜⎟ sC sC⎝⎠ =+ () () () ( )() 1232IsIsIssC sCRYs=+ 21 ⎡⎤ = ++ () () ()() () 1231XsIsRIs sCR sCR Ys=+ sC⎣⎦ 11Y => () =, pentru RC=1 (filtru trece-jos) 22Hs== 2 + 3131XsssCR sCR+ ++ () () 9 1A b) () == => A=1, amplificare la joasa frecventa Hω 222 −+ −+ ωω ω ξωω 0013 2jjω 1 0 33Q== , factorul de calitate, 1ω=pulsaţia naturala 0 2ξ 1135−± , s= c) () 1,2, 2Hs== 2 − − ()( ) 1231ss ssss++ AB Hs=+ () 12ss ss−− 1215 = = = = => () () 1AssHs=− 1ss= 535(35) 5ss− −+ −−− 12ss= 1215 = = =− =− => () () 2BssHs=− 2ss= 535(35) 5ss− −− −−+ 21ss= 125st s t−− ⎡⎤ ht e t e tσσ => () =− () () ⎣⎦ 5 xttδ= => (ytht= Răspunsul sistemului la impulsul unitate, ( )())() = => Răspunsul sistemului la () 00cosxtA tω ω=+ () ( )() {} () 00 0 0cos argyt A H t Hωω 11 =⇒= ()() 0HHωω 22 22 2 2 ωωω −+ − + ()() 0019 1 9ω 2 arg arg 1 3Hjωωω=− − + () {} {} 03ω Daca 1ω atunci () =− {} 20arg arctgHω 0 − 01ω 57jω+ 10 Pentru un SLIT răspunsul in frecventa este () Hω = 2 41jjjωω ω +++ ()() () a) Sa se determine răspunsul la impuls al sistemului h(t), b) Daţi o structura de implementare a sistemului 57 1 2ss++ a) () −+ Hs== 22 ++ 4141ssssss+ +++ () () 10 13 s+ 1 24 + + () 222Hs=− 413 13s+ ⎛⎞ ⎛⎞ ++ ss++ ⎜⎟ ⎜⎟ 24 24⎝⎠ ⎝⎠ 42 233tt t−− − σ=− + + cos 2 sinht e t e t t e t tσσ () ()()() 22 Verificati expresia raspunsului la impuls folosind mediul Matlab: >> L=-1/(s+4)+(s+2)/(s^2+s+1) >> ilaplace(L) 57 57ss++ b) () 322Hs== ++ 55441ssssss+ +++ () () b1=5, b0=7 a3=1, a2=5, a1=5, a0=4 ′′ ′ ′ sau: 5 5 4 5 7yyyyxx′′′ +++=+ 11